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Rotational structures of even-even i 48-160 ]^ nuclei are studied with the self-consistent deformed 
Hartree-Fock (HF) and angular momentum (J) projection model. Spectra of ground band, recently 
observed K — 4~, K — 5~ and a few more excited, positive and negative parity bands have been 
studied upto high spin values. Apart from these detailed electromagnetic properties (like E2, Ml 
matrix elements) of all the bands have been obtained. There is substantial agreement between 
our model calculations and available experimental data. Predictions are made about the band 
structures and electromagnetic properties of these nuclei. Some 4-qasiparticle K-isomeric bands and 
their electromagnetic properties are predicted. 

PACS numbers: 21.10.-k, 21.10.Ky, 23.20.lv, 27.70. +q 

I. INTRODUCTION 

The study of rotational structures of neutron-rich nuclei is one of the main topics of nuclear physics. For rotational 
nuclei the coupling of single-particle and collective rotational degrees of freedom leads to a variety of phenomena [1-3]. 
The neutron number N=90 marks the beginning of deformed rotational features in the rare-earth region. This is seen 
in a number of nuclei in this region [4-6]. Lighter Neodymium isotopes show vibrational like behavior, while heavier 
isotopes display more rotational- like spectra [7]. In heavier Nd isotopes the active protons are filling in sdg7/ 2 h l:1 /2 
shell above Z = 50 shell closure and active neutrons in pfhg / / 2 ii3/2 shell above N — 82 shell closure. Hence, these Nd 
nuclei are expected to show deformation properties and associated band structures. Some experimental informations 
[8-13] are now available for the band structures of neutron-rich Nd nuclei and more results can be expected with 
the existing [14-16] as well as new generation facilities [17-21] for structure studies away from the valley of stability 
towards the neutron drip line. It is thus of great interest to make microscopic study of the structure of these nuclei 
and to know quantitatively the bands to be expected for these. 

For a reliable description of the band structure, one needs a microscopic many-body self consistent procedure to give 
insight into the shape and intrinsic structure of the nucleus and angular momentum projection procedure to obtain the 
various J states from the intrinsic states. The model used should be based on residual interaction among nucleons in 
a reasonable model space. In this work we use deformed Hartree-Fock and Angular Momentum (J) projection method 
[22-24] for the study of band structures of even-even Nd nuclei (A — 148 — 160). We call the deformed HF and angular 
momentum projection the Projected Hartree-Fock (PHF) model. This self-consistent procedure ensures the correct 
particle number (for protons and neutrons) and gives eigenstates of angular momenta for the states of the bands. 
This model is known to take into account the configuration mixing in the active shell model space and hence gives 
a quantitative description of energy spectra and electromagnetic transitions in deformed nuclei [25-27]. The residual 
interaction is included self-consistently in building the deformed basis in this theory. Once the self-consistent HF 
calculation has been performed, one can build the intrinsic states of the bands by particle-hole excitations across the 
proton and neutron Fermi surfaces (the various quasiparticle configurations), besides the Hatree-Fock configuration. 
Subsequent J-projection from the deformed intrinsic configurations gives the spectra of various bands. Diagonalisation 
after projection can be done. Thus our model is very close to the shell model [28] and can predict new bands with 
higher energies and higher angular momenta. 

The neutron rich Nd nuclei are well deformed (rotational states are known in l48 Nd- 156 Nd). But the experimental 
information are quite limited. Only a few bands are known and they are not known to high spin values. No bands 
are know experimentally for the heavier neutron-rich Nd nuclei ( 158 Nd onwards). However, the current experimental 
facilities like spontaneous fission of 252 Cf, thermal neutron induced fission make it possible to study band structures of 
heavier deformed rare-earth nuclei. The experiment with spontaneous fission of 252 Cf is used by Zhang et al. [8] and 
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Gautherin et al. [10] to study yrast band, negative parity bands and bands based on K-isomers in i 52 . 154 . 156 ]^ nuclei. 
Recently, using the same technique, Simpson et al. [9] reported the two-quasiparticle (2-qp) isomeric band structure 
for 154 ' 156 Nd nuclei. They are based on 4~ isomeric state in 154 Nd and 5~ isomeric state in 156 Nd. Similar band 
structure based on 7~ isomeric state is also observed by Yeoh et al. for 152 Nd, in their recent study [29]. Motivated 
by such a good number of experimental observations and very few systematic theoretical descriptions so far of these 
nuclei we have study the band structures for a wider range of even-even Nd nuclei isotopic chain with mass numbers 
A = 148 to A = 160 in a self-consistent model (the PHF model). 

This paper is organized as follows: some details of the theoretical formalism and calculation procedure are presented 
in Section II. The results and discussion are given in Section III. The various possible intrinsic configurations used in 
our calculations are presented here. Comparisons of theoretical results with available experimental data for energy 
levels, bandhead energies and electromagnetic properties for even-even 148 - 160 Nq are made. Finally we summarize 
our work in Section IV. 



II. DEFORMED HARTREE-FOCK AND J PROJECTION METHOD 
A. Deformed HF Formalism 

The HF equation is derived from the nuclear Hamiltonian which consists of single-particle and residual two body 
interactions terms. The self-consistent HF method is a very useful technique to obtain nuclear mean held from two 
body interaction. Here we have used Surface Delta interaction (SDI) [23, 30] as two body residual interaction to get 
self-consistent wave functions. This is a reasonable interaction in this model space and has been used in many studies 
in this mass region [27, 31-34]. The residual interaction causes mixing of single-particle orbits of nucleons, leading 
to an average deformed held (Hartree-Fock field) . The deformed orbits are generated by the self consistent iterative 
of the Hartree-Fock equations. Here we use an axially symmetric basis with K quantum number for each intrinsic 
state. This is actually not a limitation of our model, because we can diagonalise among various K configurations after 
projection. 

Let rjm) be a deformed orbit which can be written as: 

= Y^Cjm\jm) (1) 

3 

where j is the angular momentum of the spherical single particle state and m its projection. The amplitude Cj m of 
the spherical state \jm) in the orbit r\ obtained by solving deformed Hartree-Fock equations 

EVC jm = e 3 C ]m + ^(hm3jim i \V\jmj 2 m 2 )p nm2M ^C] 3m3 (2) 

where ej is the single particle energy of shell model orbit j(= nlj), V denotes the two-body interaction among nucleons 
and the density matrix 

Pj2m 2 j4m 4 — (4>K \ a j i m i a 32m 2 \4>k) = ^ Cj 4 ro4^j2m 2 (3) 

rj(occupied) 

For axial symmetry one has the conditions m 3 = m and m 4 = 1712- Equations (2) and (3) are solved iteratively till the 
convergence obtained. The state \rjrn) for protons and neutrons form the deformed orbits. The Slater determinants 
for protons and neutrons constitute the deformed intrinsic state \4>k)- 

B. Angular Momentum (J) Projection 

Because of mixing in single particle orbits, intrinsic states and the multi nucleons Slater determinants are not good 
angular momentum (J) states (these are states of good K because of assumption of axial symmetries). Thus \4>k) 
does not have a unique J quantum number and is superposition of various J states: 

\4>k) =Ys C k\M ( 4 ) 
J 

To obtain the band structure and other electromagnetic properties one needs good angular momentum (J) states. 
The good-J states are obtained by J-projection from deformed intrinsic configurations. This is an important aspect 
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of our microscopic theory. Angular momentum projection from various intrinsic states gives the bands in a nucleus. 
The angular momentum projection operator is [35]: 

P k M = J dn D MK ~m R(Q) (5) 

Here R(fl) is the rotation operator e~ laJ " e~ l @ Jy e~ llJz and ft represents the Euler angles a, f3 and 7. For details of 
HF and angular momentum projection see Refs. [22, 23, 36]. The normalized state of good angular momentum (J) 
projected from intrinsic state \4>k) can be given as 

\K M ) = , P ' n4>K) (6) 
j {<t> K \P J K M \<t> K ) 

The matrix element of the nuclear Hamiltonian, containing single-particle and two-body residual interaction term, 
between two states of angular momentum J projected from intrinsic states <f> Kl and (f> K2 is: 

/ 1 1 1 r 71 

H J KlK2 = 2 (N j N j )1/2 Si n(P)d J K 1 KM(<pK 1 \He-^\ ( p K2 ) (7) 



where 



2 7+1 

N^k 2 = ^—I df3sin{l3)d J KlK2 {f}){<t>KAe-^\<l>K 2 ) (8) 



: / 7 'dp *n(J3)d J KlK2 (J3){<t>Ki\e- if)Jv \<t>K, 
J 



Euler angles a and 7 are integrated out because of axial symmetry and the remaining Euler angle (3 has to be 
integrated numerically. 

Reduced matrix elements of tensor operator T L between projected states and ipK2 are given by 



m\ TL UK.) = ^^^f^^E^ 1 / "dfi sin{P)dX{P){4>KAT^e-^ K2 ) (9) 
where the tensor operator T L denotes electromagnetic operators of multipolarity L. 



C. Bandmixing 

In general, two states \ipj£f) and \ipj^) projected from two intrinsic configurations \(t>Kx) and \<j>K 2 ) are not 
orthogonal to each other even if the intrinsic states \<j>Ki) and \4>k 2 ) are orthogonal. We orthonormalise them using 
following equation 

Y^{ H KK'-EjN J KK ,)b J K ,=Q (10) 

K' 

Here N^ K , are amplitude overlap and b J K , are the orthonormalised amplitudes, which can be identified as bandmixing 
amplitudes. The orthonormalised states are given by 

* JM = E^ M ( n ) 

K 

With these orthonormalised states we can calculate matrix elements of various tensor operators. 



III. RESULTS AND DISCUSSION 



A. The deformed single particle configurations 

The deformed HF orbits are calculated with a spherical core of 132 Sn, the model space spans the 2si/ 2 , 1^3/2, ^d 5 / 2 , 
0^7/2, O/19/2, 0h n / 2 orbits for protons and the 2p 1/2 , 2p 3/2 , I/5/2, I/7/2, 0/i 9 / 2 , 0i 13 / 2 orbits for neutrons. Single 
particle energies of these orbits used for both HF and J projection calculations are given in the Table I. We use 
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TABLE I: Single Particle Energies of the model space 



Proton 
[MeV] 


Sl/2 


d-i/2 


4/2 


97/ 2 


h-9/2 


^11/2 


3.654 


3.288 


0.731 





7.1 


2.305 


Neutron 
[MeV] 


Pl/2 


P3/2 


h/2 


iV/2 


hg/2 


*13/2 


4.462 


2.974 


3.432 





0.300 


1.487 



surface delta interaction as the residual interaction among the active nucleons in these orbits. The interaction has 
the following form (see Ref. [30]) 

V(r 12 ) = -2F(R u y 4 8 (n - R )S(r 2 - R Q )S(cosu 12 - 1) 

= -^^^'(ui^N (12) 

lm 

where R is the nuclear radius, u is the radial wave function. Here the radial wave functions are approximated to be 
same at the nuclear surface. Vq is the strength of the interaction and is taken to be 0.3 MeV for proton-proton (pp), 
neutron-neutron (nn) and proton-neutron (pn) interactions in the present calculation. This interaction gives a good 
description of the systematics of deformations in this mass region [31, 32]. 

The Hartree-Fock solutions for even-even 148 ~ 160 Nd nuclei are given in Table II. Since the prolate shape is lower 



TABLE II: Hartree-Fock solutions for 6oNd nuclei. 

A Shape E HF (Q 2 u) P ~ 

(MeV) Proton Neutron Theory Experiment [11] 



148 


Prolate 


-45 


195 


12 


275 


12 


676 





273 





201±0.004 




Oblate 


-42 


875 


-10 


745 


-10 


885 


-0 


237 






150 


Prolate 


-56 


163 


12 


332 


15 


030 





287 





285±0.002 




Oblate 


-52 


250 


-10 


962 


-12 


080 


-0 


247 






152 


Prolate 


-67 


303 


12 


387 


17 


113 





300 





349±0.012 




Oblate 


-66 


679 


-11 


139 


-14 


614 


-0 


265 






154 


Prolate 


-78 


483 


12 


600 


18 


006 





308 








Oblate 


-61 


842 


-4 


178 


-8 


337 


-0 


117 






156 


Prolate 


-91 


123 


12 


631 


20 


089 





320 








Oblate 


-82 


592 


-9 


727 


-16 


885 


-0 


256 






158 


Prolate 


-102 


957 


12 


660 


21 


614 





329 








Oblate 


-101 


640 


-11 


503 


-19 


432 


-0 


297 






160 


Prolate 


-113 


350 


12 


739 


23 


196 





338 








Oblate 


-112 


814 


-11 


574 


-19 


747 


-0 


299 







in energy than the oblate for these nuclei, we have used the prolate deformed Hartree-Fock solutions of 148 - 160 Nd for 
studying the band structures. As an example, the prolate HF orbits for 152 Nd are shown in Fig. 1. 

The quadrupole deformation parameter (/?) values obtained in the HF calculations are also listed in Table II. We 
have used the effective charges of 1.7e and 0.7e for protons and neutrons respectively. The experimental trend of f3 
values (known for A=148-152) are reproduced quite well. The j3 values show an increasing trend till A = 160. 

Because of time-reversal symmetry the nucleon orbits | ± Q) are doubly degenerate for these even-even nuclei. We 
make suitable neutrons and protons excitation across the Fermi surfaces and considered suitable intrinsic configurations 
for angular momentum projection. These intrinsic states are listed in Table III. 

As one of our main interests is to study low-lying K-isomers, so it is important to give a closer look to the levels near 
the proton and neutron Fermi surfaces. Protons orbits based on hn/ 2 (l/2) and /in/ 2 (3/2) are just below the proton 
Fermi surface and orbits based on g 7 / 2 (5/2) and /in/ 2 (5/2) are just above it. Hence, it is expected that these orbits 
can be coupled to give K = 4 and K = 5 isomeric states. In neutron side, the orbits near neutron Fermi surface are 
based on fi/ 2l hg/ 2 an d *i3/2- The K configurations from proton side is uniform but as the neutron number changes 
through out the isotopic chain there are considerable variations around the neutron Fermi surfaces leading to the 
change in the configurations of K-isomers built on quasi-neutrons. 
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FIG. 1: Prolate deformed HF orbits of 152 Nd. The orbits are denoted by ^m^. Occupied orbits are marked by (*). 



TABLE III: Intrinsic configurations for J-projection calculation for eoNd nuclei. 



Isotope Index W Configuration 



148 Bl 


0+ 


HF 




B2 


1" 


7r(3/2~ 


®5/2+) 


B3 


0" 


7r(3/2" 


® 3/2+) 


B4 


T 


77(1/2" 


®5/2+) 


B5 


4" 


vr(3/2- 


<g)5/2+) 


B6 


5" 


7r(5/2 - 


® 5/2+) 


B7 


5" 


1/(5/2" 


®5/2+) 


B8 


9+ 


7r(3/2" 


® 5/2+)® 1/(5/2^ <g) 5/2+) 


B9 


10+ 


7r(5/2" 


® 5/2+)® i/(5/2- (8) 5/2+) 


150 Bl 


0+ 


HF 




B2 


1" 


tt(3/2- 


® 5/2+) 


B3 


4" 


7r(3/2" 


® 5/2+) 


B4 


5" 


tt(5/2- 


® 5/2+) 


B5 


5" 


i/(5/2- 


<g)5/2+) 



Continued on next page. . . 
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TABLE III - Continued 



Isotope Index K" 


Configuration 








B6 


9+ 


tt(3/2- 


5/2+) 


<8)i/(5/2- 


5/2+) 




B7 


10+ 


7t(5/2" 


5/2+) 


(8)i/(5/2- 


5/2+) 


152 


Bl 


0+ 


HF 










B2 


0" 


1/(3/2+ 


(8)3/2-) 








B3 


1" 


1/(3/2+ 


(8)5/2-) 








B4 


1" 


1/(3/2" 


5/2+) 








B5 


2" 


1/(1/2+ 


(8)5/2") 








B6 


4r 


i/(3/2- 


5/2+) 








B7 


4" 


7r(3/2" 


5/2+) 








B8 


5" 


7r(5/2^ 


5/2+) 








B9 


5" 


i/(5/2- 


5/2+) 








BIO 


9+ 


7r(3/2" 


5/2+) 


®i/(5/2- 


5/2+) 




Bll 


10+ 


7t(5/2" 


5/2+) 


® i/(5/2- 


5/2+) 


154 


Bl 


0+ 


HF 










B2 


1" 


tt(3/2- 


5/2+) 








B3 


1" 


i/(3/2- 


5/2+) 








B4 


4" 


7r(3/2~ 


5/2+) 








B5 


4" 


i/(3/2- 


5/2+) 








B6 


5" 


7r(5/2~ 


5/2+) 








B7 


5" 


1/(5/2" 


5/2+) 








B8 


9+ 


7r(3/2" 


5/2+) 


®i/(5/2- 


5/2+) 




B9 


10+ 


7r(5/2" 


5/2+) 


i/(5/2- 


5/2+) 


156 


Bl 


0+ 


HF 










B2 


1" 


7r(3/2^ 


5/2+) 








B3 


4" 


7r(3/2" 


5/2+) 








B4 


5" 


7r(5/2~ 


5/2+) 








B5 


5" 


i/(5/2- 


5/2+) 








B6 


9+ 


7r(3/2^ 


5/2+) 


0i/(5/2- 


5/2+) 




B7 


10+ 


7r(5/2" 


5/2+) 


i/(5/2- 


5/2+) 


158 


Bl 


0+ 


HF 










B2 


1" 


7r(3/2~ 


5/2+) 








B3 


4" 


7r(3/2~ 


5/2+) 








B4 


5" 


7r(5/2^ 


5/2+) 








B5 


5" 


i/(5/2- 


5/2+) 








B6 


9+ 


7r(3/2" 


5/2+) 


i/(5/2" 


5/2+) 




B7 


10+ 


7r(5/2" 


5/2+) 


1/(5/2" 


5/2 T ) 


160 


Bl 


0+ 


HF 










B2 


1" 


7r(3/2^ 


5/2+) 








B3 


4" 


7r(3/2" 


5/2+) 








B4 


5" 


7r(5/2~ 


5/2+) 








B5 


4" 


1/(1/2" 


7/2+) 








B6 


8+ 


7r(3/2" 


5/2+) 


0i/(l/2- 


7/2+) 




B7 


9+ 


7r(5/2~ 


5/2+) 


0i/(l/2- 


7/2+) 



B. Energy Spectra 

Deformed Hartree-Fock and Angular Momentum Projection calculations are performed for even-even Nd isotopic 
chain with mass number A = 148 to A — 160. The energy spectra associated with each intrinsic state is obtained 
by angular momentum projection in the given model space and with surface delta residual interaction, using the 
formalism discussed in Section II. The energy spectra, band-head energies and other spectroscopic properties are 
shown in Tables IV-V and in Figures 2-14 for even-even 148 ~ 160 Nd isotopes. 



7 



TABLE IV: Comparison of calculated and experimentally observed transition energies , reduced transition probabilities 
B(E2:J J — 2), spectroscopic quadrupole moments Qs and (/-factors g(J) of 6oNd nuclei. Here B(E2) and Q(J) are calculated 
for effective charges e p — 1 + e e // and e n = e e // where e e // = 0.7. The g-factors of gf = 1.0, g\ = 0.0, <;J = 5.586 x 0.5 and 
g v a — —3.826 x 0.5 are used to calculate g(J). The experimental transition energies are taken from Ref. [7] for 148,150 Nd, Ref. 
[8] for 152 Nd and Ref. [9] for 154 ' 156 Nd. The experimental B(E2) and Qs values are from Ref. [11] and that of g(J) are from 
Ref. [12] 



Transition Transition Energy (MeV) B(E2:J — > J - 2) (e^b 2 ) Q s (eb) g(J) (nm) ~ 

J —¥ J — 2 Theory Experiment [Ref.J Theory Experiment Theory Experiment Theory Experiment 

T3S Nd 



2+ ->■ 0+ 


0.1122 


0.302 


0.5228 


0.276±0.004 


-1.4643 -1.46± 0.13 


0.6153 


0.365±0.015 


4+ ->2+ 


0.2603 


0.451 


0.7440 


0.428±0.004 


-1.8624 - 


0.6157 


0.35±0.05 


6+ ->■ 4+ 


0.4047 


0.527 


0.8136 




-2.0465 - 


0.6162 


0.267±0.050 


8+ ->■ 6+ 


0.5438 


0.576 


0.8430 





-2.1502 - 


0.6169 


_ 


10+ ->■ 8+ 


0.6765 


0.615 


0.8539 


- 


-2.2163 - 


0.6181 


-0.175±0.009 


12+ -»■ 10+ 


0.8018 


0.635 


0.8545 




-2.2599 - 


0.6192 




150 Nd 
















2+ -»■ 0+ 


0.0719 


0.130 


0.5794 


0.542±0.006 


-1.5413 -2.0± 0.5 


0.3828 


0.42±0.02 


4+ ->■ 2+ 


0.1646 


0.251 


0.8265 


0.820±0.036 


-1.9587 - 


0.3808 


0.44±0.03 


6+ -»■ 4+ 


0.2508 


0.339 


0.9074 


0.991±0.055 


-2.1506 - 


0.3777 


0.35±0.07 


8+^6+ 


0.3278 


0.409 


0.9453 


1.136±0.02 


-2.2592 - 


0.3740 


0.56±0.13 


10+ ->■ 8+ 


0.3951 


0.470 


0.9656 


0.996±0.02 


-2.3275 - 


0.3701 


0.14±0.19 


12+ -»■ 10+ 


0.4541 


0.520 


0.9754 




-2.3742 - 


0.3662 




152 Nd 
















2+ ->■ 0+ 


0.0725 


0.072 


0.6432 


0.84±0.06 


-1.6241 - 


0.3755 


_ 


4+ ->■ 2+ 


0.1683 


0.164 


0.9169 


1.04±0.05 


-2.0657 - 


0.3754 


_ 


6+^4+ 


0.2619 


0.247 


1.0069 


1.06±0.02 


-2.2705 - 


0.3755 





8+^6+ 


0.3526 


0.322 


1.0490 


_ 


-2.3876 - 


0.3755 





10+ ->■ 8+ 


0.4395 


0.389 


1.0709 


- 


-2.4619 - 


0.3755 


- 


12+ ->■ 10+ 


0.5220 


0.452 


1.0817 




-2.5124 - 


0.3756 




154 Nd 
















2+^0+ 


0.0713 


0.071 


0.6888 


0.45±0.12 


-1.6797 - 


0.3822 





4+ ->■ 2+ 


0.1633 


0.162 


0.9825 


_ 


-2.1355 - 


0.3801 





6+^4+ 


0.2485 


0.248 


1.0792 


_ 


-2.3446 - 


0.3769 


_ 


8+^6+ 


0.3243 


0.328 


1.1253 


_ 


-2.4627 - 


0.3729 


_ 


10+ -»■ 8+ 


0.3904 


0.400 


1.1499 


- 


-2.5367 - 


0.3688 


- 


12+ ->■ 10+ 


0.4482 


0.466 


1.1631 




-2.5871 - 


0.3649 




i56 Nd 
















2+ -»■ 0+ 


0.0701 


0.067 


0.7571 


_ 


-1.7615 - 


0.3740 


_ 


4+ ->■ 2+ 


0.1626 


0.155 


1.0800 


- 


-2.2409 - 


0.3739 


— 


6+^4+ 


0.2531 


0.239 


1.1863 




-2.4628 - 


0.3738 




8+ ->■ 6+ 


0.3407 


0.318 


1.2366 




-2.5893 - 


0.3737 




10+ ->■ 8+ 


0.4245 


0.391 


1.2636 




-2.6691 - 


0.3736 




12+ ->■ 10+ 


0.5043 


0.459 


1.2772 




-2.7231 - 


0.3734 




158 Nd 
















2+ ->■ 0+ 


0.0779 




0.8146 




-1.8277 - 


0.3936 




4+ ->■ 2+ 


0.1810 




1.1618 




-2.3245 - 


0.3936 




6+^4+ 


0.2823 




1.2758 




-2.5538 - 


0.3935 




8+ -> 6+ 


0.3810 




1.3295 




-2.6845 - 


0.3934 




10+ -»■ 8+ 


0.4763 




1.3578 




-2.7670 - 


0.3935 




12+ ->■ 10+ 


0.5676 




1.3722 




-2.8223 - 


0.3934 




160 Nd 
















2+^0+ 


0.0734 




0.8776 




-1.8974 - 


0.3774 




4+ ->■ 2+ 


0.1706 




1.2516 




-2.4130 - 


0.3775 




6+ ->■ 4+ 


0.2664 




1.3749 




-2.6525 - 


0.3774 





Continued on next page. . . 
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TABLE IV - Continued 



Transition Transition Energy (MeV) 


B(E2:J -»■ J - 2) (e^b* ) 


Qs (eb) 


g(J) (nm) 


J — > J — 2 Theory Experiment [Ref.J 


Theory Experiment 


Theory Experiment 


Theory Experiment 


8+ -»■ 6+ 0.3601 


1.4327 - 


-2.7893 - 


0.3775 - 


10+ -> 8+ 0.4510 - 


1.4634 - 


-2.8761 - 


0.3775 - 


12+ -> 10+ 0.5387 - 


1.4791 - 


-2.9349 - 


0.3775 - 



1. The Ground Bands 

The K = + ground bands of even-even 148_160 Nd have rotational behaviour with the J values extending upto 32ft 
in our calculations (we have shown the ground band upto about 26ft because of limitations of energy scale) with the 
excitation energy about 10 MeV. Experimentally the K = + ground bands of even-even 148 ~ 156 Nd nuclei are known 
upto about 16ft and they are quite well described by our theoretical model (see Figures 2-6). The transition energies 
of few lowest states of these nuclei in comparison with available experimental data arc compared in Table IV. For 
158 Nd and 160 Nd energy spectra are not known experimentally and we predict a good rotational band structures for 
these nuclei as shown in Figures 7 and 8 (the bands in these figures are all predictions for these two neutron-rich 
nuclei). 

2. The Negative Parity Bands 

The nucleus 148 Nd exhibits a low-lying negative parity band with 3~ state lying below 1~ in an early beta decay 
experiment by Walters et al. [37] and later it was extended by Ibbotson et al. [38]. Our band mixing (by mixing B2, 
B3 and B4 as given in Table III for 148 Nd) calculation reasonably reproduces the spectrum of this band with 3~ state 
lower in energy than 1 _ state as in experiment. The lowest band after band mixing calculation is compared with the 
experimentally observed band and shown in Fig. 2. The bandhcad energies of the calculated K = 1~ and K = 2~ 
negative parity bands are listed in Table V. 

In Fig. 3, the calculated K = 1~ band for 150 Nd is compared with the experimentally known negative parity band 
[39]. Though the calculated bandhead energy is about 400 kcV higher than that of the experimental value, good 
agreement of moment of inertia between the calculated and experimental bands help us to conclude that the negative 
parity band could be of K = 1~ (7r3/2~[541] <g> 7r5/2 + [413]) origin as predicted in Ref. [32]. We have also predicted 
the unfavoured even-J branch of this band. The intensity to populate these unfavoured states may be small, but 
still can be looked for with multi-detector arrays. The signature splitting in this band is evident from the doublet 
structure because of the rotational alignment of 7rh n / 2 . 

After mixing a few low energy intrinsic configurations (B2, B3 and B4) for 152 Nd we obtained negative parity bands. 
The calculated negative parity bands are compared with experimental K = 0~ and K = 1~ bands (see Fig. 4). The 
signature splitting is observed in our calculation for both the bands. For K = - band only lowest 3 levels of odd-J 
branch are experimentally known [40] and these are compared with our calculated favoured branch. The other branch 
is not observed in experiment. From an early experiment Hellstrom et al. [40] find & K = 2~ band but recently Yeoh 
et al. [29] assigned this as a K = 1" band with configuration j/5/2 + [642] <g> i/3/2 _ [521]. Experimentally this band is 
known upto 12ft (only even J states). We have compared our calculated K = 1~ band with this in Figure 4. 

A similar K = 1~ band is also known experimentally for 154 Nd [9]. The calculated band in comparison with available 
experimental data is shown in Fig. 5. An excellent agreement is obtained between calculated and experiment spectra 
for this isotope. No K = 1~ band has been observed experimentally for 156 458,i60j ( j c ]_ Q ur predictions for these bands 
are shown in Figures 6-8. 

3. The K-Isomer Bands 

Experimentally a low lying K = 5~ band (with a few states built on the bandhead) is known in 156 Nd [9], while for 
154 Nd only a single K = 5~ isomeric state in known [10]. We find low lying K = 5~ bands for all Nd nuclei studied in 
this article (see configurations in Table III). In our PHF calculation the K = 5^ proton 2-qp band is lower in energy 
by about 5 MeV than the K = 5~ neutron 2-qp band for 148 Nd (see Fig. 2). In all other cases (for A = 150 — 160) the 
neutron 2-qp bands are lower in energy. There is a large energy gap between the filled and next unoccupied orbit on 
the neutron side in ground state configuration of 148 Nd. Therefore it requires more energy to excite particles across 
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the neutron Fermi surface. While for the other isotopes the orbits are closed to the Fermi level so low energy neutron 
excitations are possible. There is no experimental evidence for the K-isomer band in 150 Nd. Our predictions for 2qp 
K-isomer band for 150 Nd and 152 Nd are shown in Fig. 3 and Fig. 4 respectively. The calculated K = 5~ isomeric 
band for 154 Nd is given in Fig. 5 and compared with experimental findings [10]. Our predicted K = 5~ band-head is 
about 300 keV higher than experimental value. Similar K = 5~ band is recently observed for 156 Nd [9] and this band 
is experimentally known upto 13ft. We predict reasonably this band and our calculated band is extended upto 32ft. 
Apart from these we have also studied K = AT 2-qp bands from protons/neutrons excitations and these arc listed 
in Table V. Recently, K = 4~ isomer band is observed experimentally for 154 Nd [9] with dominant configuration as 
i^5/2+[642] (g) i/3/2 _ [521]. Our PHF calculations also give a neutron 2-qp band with same configuration. Similar band 
is also obtained in recent Projected Shell Model (PSM) calculations [41]. As it can be seen from the Table V that 
proton 2-qp K = AT bands are low lying with bandhead energies less than 2 MeV (except for 152 Nd it is about 2 
MeV). For Nd isotopes with TV > 98 the i j^. (7/2) state is very close to the neutron Fermi surface so we obtained a 
low-lying neutron 2-qp K = 4 _ band with bandhead energy about 0.5 MeV for 160 Nd. 

Apart form the 2-qp bands we have also studied the 4-qp K = 8+, K = 9+ and K = 10+ bands obtained from 
the combination of proton and neutron 2-qp excitations. Our predictions for these bands are shown in Figures 2-8. 
All these bands are identical (bandhead energies and electromagnetic properties are shown in Tab. V) for the nuclei 
studied here and for neutron rich nuclei these high-K isomer bands are low-lying. The configurations for these bands 
are also listed in Table III. 

TABLE V: Comparison of calculated and experimentally observed band head energy (BHE), spectroscopic quadrupole moments 
(Qs) and magnetic moment (fx) of Nd nuclei. Here Qs are calculated with effective charges e p = 1 + e e // and e n = e e // where 
e eff = 0.7. The g-factors of = 1.0, g\ = 0.0, gl = 5.586 x 0.5 and g v 3 = -3.826 x 0.5 are used to calculate fj,. 



Nuclei K" 


BHE (MeV) 


Qs (eb) 


/i (nm) 




Theory Experiment [Ref.] 






148 Nd l-(a) s 


1.1584 0.9992 [38] 


-1.2078 


1.7655 


2" (a) 


1.8100 


1.2936 


2.0305 


4- (a) 


1.6189 


2.4688 


3.6675 


5- (a) 


3.2242 


2.6836 


4.7003 


5"(b) 


8.3077 


2.7182 


0.2112 


9+ 


9.1619 


3.2392 


3.8401 


10+ 


10.2181 - 


3.2406 


4.8406 



l-(a) 


1 


.2541 


0.8529 [39] 


0.5138 


0.3045 


4- (a) 


1. 


.6754 




2.6097 


3.5146 


5"(b) 


2 


.6951 




3.0178 


0.2163 


5- (a) 


3 


.1465 




2.8669 


4.5417 


9+ 


4 


.1121 




3.6184 


3.8383 


10+ 


5 


.1456 




3.6343 


4.8607 


0~(b) § 


1 


.1947 


1.1486 [40] 


-0.9279 


0.3723 


l-(b) 


1. 


.4778 




0.5629 


0.2091 


4- (a) 


2 


.0516 




2.7567 


3.4989 


4"(b) 


1. 


.1094 




2.6080 


0.1386 


5"(b) 


2 


.3347 




2.2677 


0.1347 


5- (a) 


3 


.7741 




3.0344 


4.5384 


9+ 


3 


.9071 




3.8367 


3.7329 


10+ 


5 


.0516 




3.8698 


4.7339 


l-(b) 





.9412 


0.9619 [9] 


0.5625 


0.2993 


4- (a) 


2 


.8622 




2.8575 


3.5161 


4"(b) 


1. 


.0357 


1.2980 [9] 


3.0174 


0.2687 


5"(b) 


1. 


.6644 


1.348 [10] 


3.3648 


0.1306 


5- (a) 


3 


.1832 




3.1466 


4.5386 


9+ 


3 


.0651 




4.0577 


3.7488 


10+ 


4 


.0333 




4.0969 


4.7390 



1" 


(a) 


1.8127 


3.0021 


3.5005 


4" 


(a) 


1.7308 


3.0021 


3.5005 


5" 


(b) 


1.2186 1.4312 [9] 


3.5664 


0.2694 



Continued on next page. . . 
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TABLE V - Continued 



Nuclei 




BHE (MeV) 
Theory Experiment |Ref.] 


Qs (eb) 


/! (nm) 




5" (a) 


3.1181 


3.3117 


4.5353 




9+ 


3.1211 


4.3061 


3.8841 




in+ 


A 409^ — 


A ^990 




158 Nd 


l-(a) 


2.0440 


0.6138 


0.2937 




4- (a) 


1.9652 


3.1179 


3.5127 




5"(b) 


2.8203 


3.6161 


0.2186 




5- (a) 


3.3823 


3.4399 


4.5418 




9+ 


5.5992 


4.4629 


3.8729 




10 


O.0933 - 


A A AO*~r 

4. 4437 


4.8306 


160 Nd 


l-(a) 


1.7625 


0.6379 


0.3009 




4"(b) 


0.5272 


3.3536 


0.2544 




4- (a) 


1.6855 


3.2417 


3.5023 




5- (a) 


2.9253 


3.5807 


4.5335 




8+ 


2.2104 


4.4304 


3.8387 




9+ 


3.2361 


4.4830 


4.8425 



(a) proton excitation (b) neutron excitation 
$ lowest state is 3~ 

§ K = 0~ band, Qs and [i are given for 1~ state 



C. Electromagnetic Properties of Bands 

1. B(E2) 

Apart from the energy spectra we have calculated the reduced transition matrix elements e.g. B(E2), B(M1). The 
B(E2) values for a transition from an initial state a.J\ to final state /3J2 is given by 



B(E2; aJi -> j3J 2 ) 



1 



2Ji + l 



i—p,n 



mi 1 



(13) 



where the summation is for quadrupole moment operators of active protons and neutrons. The effective charges of 
1.7e and 0.7e are used for protons and neutrons respectively in our calculations. The calculated B(E2) values for 
ground state bands are given in Table IV. Also, in Figures 9 and 10 we have shown the behavior of B(E2) values 
for ground bands, negative parity bands as well as K-isomer bands for the nuclei studied in this paper. From these 
figures we find that for the ground band of all the Nd isotopes considered, B(E2) values are smoothly varying with 
increasing J. The B(E2; J — > J — 2) values for the ground band show similar trend through out the isotopic chain 
studied here. It is maximum in the range J=12 to J=16 and falling down beyond this range. For other bands in case of 
J — > J — 2 transitions the values are increasing with J and after a certain spin its nearly constant. For the transitions 
J — > J — 1, the maximum B(E2) values are around 1.0 (eb) 2 for 2-qp bands and a little bit lower for 4-qp bands. 
The B(E2) values for K = 1~ bands(proton excitation) are calculated and displayed in Fig. 11. Similar variations 
in B(E2) values are observed for all the nuclei considered in our calculations. The B{E2;2 + — > + ) systematic is 
compared with experimentally available values. In Figure 12, the experimental trend reasonably reproduced with our 
calculations. 

2. B(M1) 

The reduced magnetic dipole transition moments between initial and final states are given by 



B(M\;aJ x /3 J 2 ) = 



1 3 

2 Ji + 1 4^ 



T,^K:\\9ik+gi Sl \\n J 1 i ) 



(14) 



11 



where gi and g s are orbital and spin g-factor respectively. The g-factors of gi = and g s — 5.586 x 0.5^at for 

protons and gi — O^n and g s = —3.826 x 0.5/xjy for neutrons are used for B(M1) calculations. The quenching of spin 
g-factors by 0.5 is taken in account to consider the core polarization effect [42, 43]. 

B(M1) values for K-isomer bands are shown in Fig 13. For K = 5~ proton 2-qp band B(M1) values are gradually 
increasing with J for all nuclei. B(M1) values for neutron 2-qp band in 14 8, 150,152,154^ also show a same behaviour. 
But for 156 ,i58^^ ^ B(M1) values are constant at higher J. K = AT neutron 2-qp band in 160 Nd also shows a 
similar trend. From Fig. 13, for 148 Nd we can see that at high spin above J=25, both the 2-qp bands show staggering 
in B(M1) values and it is more prominent in case of neutron 2-qp band because of the doublet structure of these 
bands at high spin values. For transitions within doublets, the terms contributing to reduce Ml matrix elements 
are large resulting this staggering effect. For 4-qp bands we notice that the B(M1) values are very small, where the 
contributions from neutron spin is comparable with that of proton; but as g-factor of neutrons is opposite sign to 
that of protons, they cancel each other giving a very small net quantity. Similarly for K = 1~ bands (from proton 
excitations) the B(M1) values are also very small and shown in Fig. 14. 



3. Spectroscopic Quadrupole Moment (Qs) 
The spectroscopic quadrupole moment of a state with angular momentum J is given by 

Qs{J) = (2J+ 1)V2 (T ) 1/2 ^^H £ <&W*k) ( 15 ) 

We have calculated the electric quadrupole moments for band head as well as other spin states for all bands. The 
calculated values are given in Tables IV-V. The experimental values are only known for ground state band and 
comparisons are given in Table IV. The known experimental results are nicely reproduced in our calculations. 



4- Magnetic Dipole Moment (fi) and g-factor 

The magnetic dipole moment (//) of a state with angular momentum J can be expressed as 

MJ) = (2jTW CJj ° Jj I ^ ^KMh + glsiWK)) (16) 



where i=p and n for protons and neutrons respectively. 
The gyromagnetic factor (g-factor) of state J is defined as 

9(J) - ^ (17) 

where fi(J) is the magnetic moment of state J. 

The g-factors for ground state band for different J values are calculated and are given in Table IV. Experimental 
data are only available for 148 > 150 Nd [11]. In our calculations, for 148 Nd the values are overestimated and for 150 Nd, the 
calculated values are reasonable in comparison to experimental values except for J=10. The magnetic dipole moments 
for K-isomer and negative parity bands are calculated and given in Table V. No experimental data is available for 
these bands. The magnetic moments for the K = 5~ neutron 2-qp bands are lower than the proton 2-qp bands. 



IV. CONCLUSION 



The neutron-rich even mass Nd nuclei have interesting quadrupole deformation properties and band structures. We 
have presented results for such nuclei, using the microscopic model of angular momentum projected deformed Hartrec- 
Fock method. Some of our results correlate well with known experimental findings about the deformation and the 
spectra of these nuclei. Using our PHF model with band mixing we predict many bands (including K-isomer bands) 
upto high spin values. The spectra for both positive and negative parity, two-quasiparticle and four-quasiparticle 
bands as well as electromagnetic properties like B(E2), B(M1) of these bands are obtained. The known systematics 
of the band structures are given quite well with this model. These results can be useful in future experimental studies 
of these nuclei. As a whole our microscopic model can be very useful to study the band structure of nuclei and 
electromagnetic properties in this mass region to high spin values. 
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FIG. 2: Comparison of theoretical spectra with experimental spectrum of various band for 148 Nd. (a) Ref. [7], (b) 

Ref. [38]. 
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FIG. 3: Same as Fig. 2 but for 150 Nd. (a) Ref. [7], (b) Ref. [39]. The K = 1" band of our calculation has both 
odd- J and even- J (unfavoured) states (see discussion in Section III B). 
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FIG. 4: Same as Fig. 2 but for 152 Nd. (a) Ref. [8], (b) Ref. [40]. 
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FIG. 5: Same as Fig. 2 but for 154 Nd. (a) Ref. [9], (b) Ref. [10]. 
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FIG. 6: Same as Fig. 2 but for 156 Nd. (a) Ref. [9]. 
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FIG. 7: Theoretical energy spectra of 158 Nd. 
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FIG. 8: Theoretical energy spectra of 160 Nd. 
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FIG. 9: B(E2) values for various bands of 6oNd nuclei. 



18 




FIG. 10: B(E2) values for various bands of 6oNd nuclei. 
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FIG. 11: B(E2) values for K = 1" band(proton 
excitations) of 6oNd nuclei. 
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FIG. 12: Variation of B(E2;2+ ->• 0+) for 60 Nd 
nuclei. The experimental values are taken from 
Ref. [11]. 




FIG. 13: B(M1) values for 60 Nd nuclei. 
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FIG. 14: B(M1) values for K = 1 band(proton excitations) of 6oNd nuclei. 



